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3 Abstract. We seek to reconcile observations of small source sizes in the 

4 solar corona at 327 MHz with predictions of scattering models that incor- 
s porate refractive index effects, inner scale effects and a spherically diverg- 
e ing wavefront. We use an empirical prescription for the turbulence ampli- 

? tude Cjj(R) based on VLB1 observations by Spangler and coworkers of com- 

a pact radio sources against the solar wind for heliocentric distances R ~ 

o 10-50 R Q . We use the Coles & Harmon model for the inner scale k(R), that 

10 is presumed to arise from cyclotron damping. In view of the prevalent un- 

11 certainty in the power law index that characterizes solar wind turbulence at 

12 various heliocentric distances, we retain this index as a free parameter. We 

o find that the inclusion of spherical divergence effects suppresses the predicted 
14 source size substantially. We also find that inner scale effects significantly 
is reduce the predicted source size. An important general finding for solar sources 
i6 is that the calculations substantially underpredict the observed source size. 
I? Three possible, non-exclusive, interpretations of this general result are pro- 
is posed. First and simplest, future observations with better angular resolution 
w will detect much smaller sources. Consistent with this, previous observations 

20 of small sources in the corona at metric wavelengths are limited by the in- 

21 strument resolution. Second, the spatially- varying level of turbulence C%(R) 

22 is much larger in the inner corona than predicted by straightforward extrap- 

23 olation Sunwards of the empirical prescription, which was based on obser- 

24 vations between 10-50 R Q . Either the functional form or the constant of pro- 
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25 portionality could be different. Third, perhaps the inner scale is smaller than 

26 the model, leading to increased scattering. 
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1. Introduction 

27 Refractive scattering of radiation by density turbulence in the Sun's corona and so- 

28 lar wind leads to angular broadening of embedded radio sources, and of cosmic sources 

29 observed through these media. The process is similar to the twinkling of stars and mod- 

30 ified "seeing" caused by density turbulence in Earth's atmosphere and ionosphere. This 

31 scattering process has been investigated for many years using geometrical optics [e.g., 

32 Steinberg et al 1971] and the parabolic wave equation [e.g., Lee and Jokipii, 1975; Coles 

33 and Harmon 1989; Bastian 1994; Cairns 1998]. 

34 Scattering is thought to affect the observed properties of type II and III solar radio 

35 bursts in several ways: greatly increasing the angular sizes of the sources [e.g., Riddle 
se 1974], causing the time profiles to have exponential decreases [e.g., Robinson and Cairns, 

37 1998], and causing anomalously low brightness temperatures at decametric wavelengths 

38 [e.g., Thejappa & Kundu 1992]. 

39 The primary motivation of this paper is to investigate the constraints imposed on models 
* of density turbulence in the solar corona by recent observations at 327 MHz [Mercier et al., 

41 2006], made by combining visibilities from the Giant Metrewave Radio Telescope (GMRT) 

42 in Pune, India, and the Nancay Radioheliograph (NRH) in France. The maps of Mercier 

43 et al. [2006] show structures ranging from the smallest observed size of 49" to that of the 

44 whole Sun, with dynamic ranges as high as a few hundred. These features make them the 

45 best meter wavelength snapshot maps of the solar corona to date. Mercier et al. [2006] 

46 found the smallest steady angular size of type I solar noise storms to be 49" in their high 

47 dynamic range, full disk, 17 second snapshots. We therefore adopt the smallest observed 
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source size of 49" in these maps to be a canonical number for comparison with our model 
predictions. The only other two-dimensional map showing small source sizes that we are 
aware of is that of Zlobec et al [1992], who observed a source as small as 30 at 327 MHz. 
However, the dynamic range of their map was severely limited, and included only a very 
small range of size scales. It should be noted that the lower limit to the observed size is 
imposed by the resolution of the instrument if scattering by density turbulence is weak 
enough. The scattering calculations shown below imply that there is a possibility that 
smaller solar sources will be detected in the future by instruments with improved angular 
resolution. 

In this paper we use a formalism based on the paraxial wave equation and the structure 
function, together with observationally based models for the density turbulence that scat- 
ter the radiation, to predict the size of sources in the solar corona at 327 MHz. Our results 
can be interpreted as the scatter-broadened image of an ideal point source in the solar 
corona. We have used an empirical model for the amplitude C%{R) of coronal turbulence 
that is based directly a fit to the scattering measure obtained from VLBI observations of 
cosmic sources broadened by scattering in the outer solar corona and inner solar wind. 
Here R is the heliocentric distance. We have assumed that this model is valid through- 
out the corona, specifically at smaller R. We also consider the effects of spherical and 
plane wave propagation, variations of the inner scale k(R) and power-law index a of the 
turbulence on the predicted source sizes. In most cases, we find that the models predict 
sizes that are at least an order of magnitude below the smallest observed size of 49 at 
327 MHz. Our formalism and analyses differ primarily from those of Bastian [1994] in the 
models for C%(R) and the electron density profile n e (R), while our applications are to 



DRAFT 



February 8, 2011, 4:13am 



DRAFT 



X - 6 SUBRAMANIAN & CAIRNS: CONSTRAINING CORONAL TURBULENCE WITH RADIO DATA 

?! metric rather than centimetric and decimetric emissions. Since our predictions are much 

72 smaller and Bastian's [1994] predictions much larger than 49" at 327 MHz, the analyses 

73 demonstrate the importance of knowing C%(R), n e (R), and k(R) much better for future 

74 observations and predictions of solar sources. These quantities are also relevant to the 

75 heating and outward flow of the coronal plasma, with activity localized to specific ranges 

76 of R potentially leading to larger Cjj(R) and so enhanced scattering at, say, decimetric 

77 frequencies than expected at, say, metric frequencies. 

78 The paper is organized as follows. In § 2 we summarize the scattering formalism and 

79 observations of the density turbulence. Coronal density models are described in § 3. The 
so results are presented in § 4, including estimates of the predicted angular broadening, and 
8i the implications for coronal density turbulence. The conclusions are presented in § 5. 

2. Angular broadening 

We first consider the angular broadening predicted by the empirical formula of Erickson 



82 Here A is the observing (free-space) wavelength in meters and D is the elongation in 
ss units of R Q . If we take D = 1.056i?©, which is where 327 MHz emission would originate 



ss much larger (60 times) than the observed 49 . It points to a significant difference between 
se the situation for observations of celestial background sources against the solar wind, for 
s? which Erickson's (1964) formula is well accepted, and observations of solar radio events 



(1964): 




arcminutes . 



(1) 



84 



according to the hybrid density model described below, then 9 ~ 50' for 327 MHz, which is 
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as that originate in the solar corona. This difference will be specifically addressed in the 
as Discussion section below. 

2.1. Density turbulence 

Density turbulence in the Sun's corona and solar wind is modeled here by writing the 
three-dimensional isotropic spatial power spectrum S n (k, R) of the fluctuating part 5n of 
the electron density n e as [cf., Lee & Jokipii 1975; Rickett 1977; Coles and Harmon, 1989; 
Bastian, 1994; Cairns 1998; Spangler, 2002] 

S n (k, R) = ((5n) 2 }(k,R) = C 2 N {R)k~ a e~ k2 ^ . (2) 

to Here k and R are the (isotropic) wavenumber and radial distance (in units of Rq), re- 

91 spectively, C%(R) models the level of turbulence, a is the power-law index, and is the 

92 wavenumber corresponding to the inner scale of the turbulence. While it is fairly well 

93 established that the turbulence spectrum largely follows the Kolmogorov scaling (with 

94 a = 11/3) at scales larger than about 100 km, there is some evidence that it flattens, 

95 with a decreasing to values as low as 3, at scales between a few km and a few hundred 

96 km [Bastian, 1994]. There is also some evidence for variation of the turbulence power 

97 law spectrum with heliocentric distance [e.g., Efimov et al., 2008]. Furthermore, there 

98 is evidence for significant variation in the index between the slow and fast solar wind 

99 [Manoharan et al., 1994]. We therefore retain a as a parameter. It may be noted that 

100 some authors use a power law index of 5/3 to describe the one-dimensional Kolmogorov 

101 spectrum; the index they refer to is equal to a — 2. 

The empirical model we use for C%(R) was originally mooted by Armstrong & Woo 
[1980] and later refined, based on VLBI observations between 10-50 R Q , by Spangler & 
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Sakurai [1995] and Spangler [2002] among others: 



C 2 N (R) = 1.8 x 10 



10 



(■ 



R 



-3.66 



(3) 



\10R Q 

102 The dimensions of C% depend on a, being m~ a_3 . The normalizations for C% differ by 

103 about a factor of a few and the power-law index with R ranges from —3.66 to —4 in these 

104 works, presumably due to solar wind variability. 

The inner scale U is modeled using Coles & Harmon's [1989] model which agrees roughly 
with their observations, 



qi (R) = ni(R)/3V A (R) = 2n/k(R) 



2ir 



684 ne(R)- 1 / 2 



knT 1 



(4) 



los where is the ion cyclotron frequency, Va is the Alfven speed and n e is the electron 

106 density in cm -3 . This model is interpreted conventionally in terms of cyclotron damping 

107 by MHD waves. We use this definition for the inner scale throughout this paper, except 
los in two cases where we artificially set U equal to a very small value. 

mo A popular alternative prescription for C% supposes that C% is oc the square of the 

no background electron density. Such a prescription has a constant of proportionality, which 

in is often determined via observed values of the phase structure function (e.g., Bastian 

U2 1994). The magnitude of the phase structure function in turn, is very dependent on the 

us elongation to which it is referenced. We discuss this issue further in § 4. 

2.2. Plane vs spherical wave propagation 

n4 Scattering depends quantitatively on whether the wavefront is planar (1-D) or spherical 

us (3-D). When a source is embedded in the scattering medium, it is often appropriate to 

no adopt a formalism that includes the spherically diverging nature of the wavefront. The 

a? geometry for spherically diverging propagation is shown in Figure 1. 



Figure 1 
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Similarly, when a plane wave illuminates the scattering medium, a 1-D planar formalism 
is standard. In this case an observer is typically sensitive only to scattering regions (eddies) 
with sizes of order the baseline length s. In the spherically diverging situation, however, the 
observer is sensitive to a range of eddy sizes given by sa/b, where a is the (continuously 
varying) distance of the scattering screen from the source and b is the distance of the 
observer from the source. In our situation, this is tantamount to saying that the effective 
baseline for spherical wave propagation is [Ishimaru, 1978] 

s cff = sR/{R 1 - Ro) . (5) 

us This is the basic difference between Eqs (8) and (9) discussed below. 

no In the solar situation, radiation from an embedded coronal source is subject to scat- 

120 tering as it propagates to the observer. Since the radiation is generated near f p and 2f p , 

121 scattering effects are expected to be largest in the source and in its vicinity, assuming 

122 that f p (R) decreases monotonically with increasing R. Spherical effects are expected to 

123 arise in two ways. Firstly, scattering will maximally distort an initially plane wavefront 

124 close to and in the source. Secondly, on a larger scale, the solar wind density is expected 

125 to be spherically symmetric, with radiation being refracted towards the radial direction. 

126 Accordingly, spherical divergence effects are expected to be vital. They are explicitly 

127 calculated below and shown to be quantitatively important. In contrast, the planar for- 

128 malism is expected to be appropriate when the source, scattering region(s), and observer 

129 are all far apart, as assumed in calculations for pulsars and other celestial sources. 
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2.3. Structure function 

The starting points for the expression we use for the scattering angle are equations 

(4)-(7) of Coles et al [1987] that specify the structure function and the mutual coherence 

function using the parabolic wave equation (PWE) formalism that includes small-angle 

refractive scattering and diffraction, but not reflection. For the sake of completeness, 

we reproduce them below. The asymptotic forms of the gradient of the phase structure 

function D(s,R) are 



M D{s^R) = ^ 2 T{l-^^)cl{R)rl\H l {R) a ~'s\ for s cS « k(R) , (7) 



where r e is the classical electron radius, A is the observing wavelength and s c g is the 
effective interferometer spacing. It is noted that this formalism is valid only for 2 < a < 4; 
in particular, equations (6) and (7) diverge at a = 4 owing to the behavior of the term 
T(l — (a — 2)/2). It may also be noted that the branches (6) and (7) do not meet at 
s cff = k; the ratio of (7) to (6) is equal to (l/2)(a - 2) {li/s) a ~ A Y{\ + (a - 2)/2), and at 
Scff = h this is equal to unity only for a = 4. 

The effective interferometer spacing s cff is equal to s for the case of plane wave propa- 
gation, but is equal to sR/(R\ — Rq) for spherical wave propagation, as discussed in § 2.2 
and Figure 1. The phase structure function for the cases of plane wave and spherical wave 
propagation are 



d_ 
M 



D(s,R) 



8tt 2 r(l-(a-2)/2) 



C 2 N {R) t\ A 2 s a ~ 2 , fors cfr >/,(/?), (6) 



2"- 2 (a -2) r(l + (a-2)/2) 




for plane wave propagation , 



(8) 




for spherical wave propagation , 



(9) 
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136 where the lower limit of integration Ro is the radial distance from which scattering is 

137 assumed to be effective (we take this to be equal to the fundamental emission level), and 
us the upper limit Ri corresponds to the observer (here at R± — 1 AU). All quantities are 
139 assumed to have spherical symmetry and the path is assumed to be radial. 

Scattering depends sensitively on the ratio of the radiation frequency / to the local 
electron plasma frequency f p (R) [Cairns, 1998]. Equations (16) and (22) of Cairns [1998] 
include the effects on refractive scattering that arise from f p (R) being non-zero and varying 
with position between the source and observer. By analogy with these equations we write 

D » As} = £ i-mr-/r- m D(s - R)dR (10) 

for plane wave propagation and 

wo for spherical wave propagation, respectively. 

The scattering angle is conventionally defined using a coherence scale s in the following 
manner [e.g., Coles et al., 1987; Bastian, 1994]: 

9 C = (2ns /\)~ l , (12) 

where 

D,(s ) = 1 (13) 

i4i and D*(s) is either equal to D p f(s), defined by (10), or D s f(s), defined by (11), in ap- 
«2 propriate limits. This scattering angle 9 C can be interpreted as the predicted size of an 
143 idealized point source. 
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2.4. Density Models 

A model for n e (R) in the corona and solar wind is required to be able to predict the 
angular broadening. The density model is required for computing the inner scale, which is 
defined in the next subsection. Since there is no universally accepted model, we initially 
consider four representative density models. One is the four-fold Newkirk density model 
for the corona, based on eclipse observations [Newkirk, 1961]: 

n 4n (R) = 4 x 4.2 x 10 4 x 10 4 - 32/i? cm" 3 , (14) 

The second model is derived from the frequency drift rate of interplanetary type III bursts 
[Leblanc et al., 1998]: 

n lh (R) = 3.3 x 10 5 iT 2 + 4.1 x 10 6 iT 4 + 8 x 10 7 iT 6 cm" 3 . (15) 

The third model considered is due to Aschwanden et al. [1995]. It is based on the drift 
rates of type III bursts [Alvarez & Haddock, 1973] in the outer corona and solar wind 
(/ < 10 MHz) and assumes an isothermal barometric atmosphere for the lower corona: 

mfe) \ R> l + R 2 



(16) 

n Q exp I - 2=± ) , R < 1 + R 2 



n a (R) = < 

where p = 2.38, Uq = 4.6 x 10 8 cm -3 , n\ = Uq exp(— p), p = 0.1 and R 2 = pp.. The fourth 
model is a "hybrid", using the Aschwanden & Benz [1995] model for the lower corona and 
the four-fold Newkirk model multiplied by a normalization factor (to ensure continuity) 
in the upper corona. In other words, the density nh y h(R) of the hybrid model is 

n ^ h{R) -\n a (R), R<l + R 2) (17) 
144 where A = 0.324 is the normalization factor that ensures continuity. 
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145 Figure 2 shows f p (R) = 8.97n e (i?) 1 / 2 kHz for all four density models, with n e in units 

we of cm~ 3 . 

147 The figure shows that the highest frequency predicted for R > 1 by the Leblanc et al. 

148 [1998] density model (Eq 15) is less than 100 MHz. Since our observing frequency is 327 

149 MHz, this model is therefore unsuitable for our purposes. However, the four-fold Newkirk 

150 model (14), Aschwanden & Benz [1995] model (16), and hybrid model can account for 

151 f P — 327 MHz for R > 1. However, since the Aschwanden & Benz [1995] model predicts 

152 unrealistically low densities (and consequently f p ) for R > 1 + R 2 — 1.23, only the hybrid 

153 model is considered further below. Fundamental emission at 327 MHz emanates from 

154 a heliocentric distance of 1.055 R Q with this model. In order to avoid the singularity 

155 in the integrand in Eqs (10) and (11), we start the integration at R = 1.056 R Q . In 

156 other words, we start the integration from a distance of approximately 700 km above 

157 the height at which 327 MHz fundamental emission originates. This distance is smaller 

158 than that corresponding to the frequency difference Af between the minimum frequency 

159 of fundamental emission at a given location and the local value of f p , so that avoiding 
mo the singularity is correct. This positive frequency difference exists because conservation 

161 of energy and the wave dispersion relations force the standard nonlinear Langmuir wave 

162 processes for fundamental and harmonic emission to produce radiation of order several 

163 percent above f p and 2f p [e.g., Cairns, 1987a,b], with the value of A/// p depending on 

164 the beam and plasma parameters. 

2.5. Inner scale effects 

165 We next discuss the need for including inner scale effects in our treatment. In general, 

166 inner scale effects are important if the baseline is smaller than the inner scale. 
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Figure 3 shows the inner scale (in km) given by Eq (4), as a function of heliocentric 
distance for some of the density models discussed in the preceding section. Clearly, the 
inner scale is quite dependent upon the density model. The inner scale for the Aschwanden 
& Benz [1995] model far exceeds those for the other models, since the density with this 
model for R > 1.23 in unrealistically low; we have therefore chosen not to depict the k 
for this model in Figure 3. As explained in the preceding section, we choose to use only 
the hybrid model from now on, for it is by far the most realistic one. 

In order to ascertain the importance of inner scale effects, we compare k with the longest 
baseline (that determines the smallest source size), assuming that the source is situated 
at the fundamental emission level for 327 MHz, and the observer is at 1 AU. In order 
to compute the longest baseline, we set 49" = 1.22A/s, where 49" is the observed source 
size and A is the (free space) observing wavelength (1 meter). This yields an effective 
baseline of s ~ 5 km. In order to ascertain the relevance of the inner scale (i.e., whether 
we should be using Eq 6 or 7) we need to compare the effective longest baseline s e ff with 
the inner scale. It may be noted that a typical interferometer measurement involves a 
range of baselines, and while the longest baseline we have computed above is the one that 
limits the smallest observable source size, baselines shorter than this one do contribute 
to the overall measurement. However, our approach is appropriate because the longest 
baseline for a given source size is the largest length scale at which there is appreciable 
power. Inner scale effects are relevant only for the branch for which the baseline is <C the 
inner scale, as in Eq 7. If the longest relevant baseline is smaller than the inner scale, then 
it follows that the rest of the baselines in the problem automatically satisfy this criterion. 
Our approach thus provides a useful estimate of the importance of inner scale effects. 



Figure 3 
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wo As discussed above, s e ff = s for plane wave propagation, while s c g = sR/{R\ — Rq) for 

loi spherical wave propagation. We show the ratio of s cS to Zj in Figure 4. 

192 When considering plane wave propagation, s cS = s, and the solid line in Figure 4 

193 shows that s e ff is mostly > Zj. For heliocentric distances greater than about 30 R Q , s e ff 

194 does become somewhat smaller than Zj, but most of the contribution to the scattering 

19 5 kernel arises from distances well inside 30 R Q . We should therefore use Eqs (6) and 
«6 (10) for plane wave propagation. On the other hand, when considering spherical wave 
197 propagation, s cff = sR/{R\ — Rq), and the dotted line in Figure 4 shows that s c g is < U 
us for all R. The appropriate equations to use for spherical wave propagation are therefore 

199 Eqs (7) and (11). 

3. Results 

3.1. Plane wave propagation 

200 We first consider plane wave propagation, which is more appropriate for waves emanat- 

201 ing from a background object that is far from the scattering medium. 

202 3.1.1. Sgff ^> Zj 

203 Since Figure 4 shows that s c g > Zj for plane wave propagation for R < 30i? Q , the 

204 appropriate branch to use is Eq (6). At heliocentric distances greater than about 30 R Q , 

205 s c fj becomes marginally less than the inner scale, but we have verified numerically that 

206 this is immaterial, since most of the contribution to 9 C takes place well within 30 R Q . 

207 Figure 5 uses Eqs (6), (10) and (13) to predict the scattering angle 9 C (1 AU) at the 

208 Earth, in arcseconds, as a function of the power law index a for plane wave propagation 

209 at / = 327 MHz. The thin line is for fundamental emission and the thick line is for second 

210 harmonic emission. Removal of the refractive index effect in Eq (10), meaning the factor 
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a of (1 — fp/ f 2 )^ 1 ), causes a negligible change in the result. The predicted source size is 

12 slightly smaller for second harmonic emission, 

is 3.1.2. s eff < U 

14 Although Figure 4 demonstrates that s c s > U for plane wave propagation, we never- 

15 theless investigate the predicted scattering angle for plane wave propagation, while using 

16 branch (7). The results are shown in Figure 6. 

17 Evidently, the source sizes predicted approach the observed size for relatively steep 
is turbulence spectra; for spectra that are steeper than Kolmogorov (i.e., a > 11/3), the 
is predicted source sizes exceed the observed one. In order to investigate inner scale effects 

20 we compute the scattering angle for plane wave propagation with the inner scale set to 

21 an artificially low value of 1 m, instead of being computed self-consistently from Eq 4. 

22 The results are shown using the heavy lines in Figure 7. When inner scale effects are thus 

23 removed, it is clear that the predicted source size increases, especially for flatter spectra. 

3.2. Spherical wave propagation 

24 As discussed earlier, spherical wave propagation is appropriate when the source is em- 

25 bedded in the scattering medium, as is the case here. 

26 3.2.1. S e ff li 

27 Since Figure 4 shows that s cS <C l { for spherical wave propagation, the appropriate 
23 branch to use is (7). 

29 The solid line in Figure 7 predicts 6 C {\ AU) for spherical wave propagation, using (7), 

30 (10) and (13). Clearly, the predicted scattering angle is at least 25 times smaller than 

31 the observed one. The dashed line, on the other hand, is computed by artificially setting 

32 U = lm, while still using (7). This is tantamount to neglecting inner scale effects. We 
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233 observe that for flat spectra, inner scale effects substantially reduce (by over an order 

234 of magnitude) the predicted 9 C , but that this difference is progressively reduced as a 

235 increases. There are negligible differences between the results for fundamental and second 

236 harmonic emission, and removal of the refractive index effect causes a negligible change 

237 tOO. 

238 3.2.2. S e ff ^> li 

239 In keeping with the spirit of our treatment for plane wave propagation, we investigate 
2« the predicted scattering angle for spherical wave propagation while using branch (6), 

241 which assumes that s e ff > U. We do this despite Figure 4's prediction that s e ff is < k for 

242 spherical wave propagation. 

243 The difference between fundamental and second harmonic emission are negligible. 

244 3.2.3. Spherical vs plane wave propagation 

245 Although we have investigated several different cases, our attention has been focussed 

246 mainly on two issues: first, the difference between the source sizes predicted for plane 

247 wave and spherical wave propagation, and, second, the influence of the inner scale. We 

248 now compare the plane wave and spherical wave results directly in Figure 9, assuming 

249 fundamental emission, and employing an inner scale that is computed self-consistently 

250 using Eq (4). 

251 It is clearly evident from Figure 9 that spherical divergence effects decrease the predicted 

252 scattering angle by around two orders of magnitude as compared to the plane wave case. 

253 This is best seen by comparing the same branch (say, s e s > U) for the plane wave and 

254 spherical wave cases. Thus, spherical divergence effects should be quantitatively important 

255 for scattered solar radio emission and plane wave results should be used with great caution. 
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4. Discussion and Summary 

256 The highest resolution meter wavelength observations of the solar corona reveal com- 

257 pact sources around 49" in size at 327 MHz. The main aim of this paper is to employ 

258 an observationally-motivated model for the turbulence amplitude C^r, and see what it 

259 implies for the predicted scattering angle for radio sources located in the solar corona. 

260 We reference our calculations to the same frequency (viz. 327 MHz) at which the small- 

261 est source size is observed. We employ the parabolic wave equation, together with the 

262 standard asymptotic forms for the phase structure function, which are valid for situations 

263 where the effective baseline is either much larger or much smaller than the inner scale. 

264 We define the predicted scattering angle 6 C via Eq. (12) as the angle where the phase 

265 structure function falls to 1/e times its peak value. Effectively, this means that the scat- 

266 tering angles predicted here should be interpreted as the scatter-broadened image of an 
26? ideal point source in the solar corona. The real source will have an intrinsic size (i.e., it 

268 will not be a point source) and the observable source will be the convolution of the in- 

269 trinsic source profile with 9 C , provided there are no instrumental limitations. The results 

270 in this paper should therefore be regarded as lower limits to the observable source size 

271 set by scattering. The general consensus now seems to be that there is not much about 

272 the instrinsic source size that can be gleaned from scatter-broadened images [Bougeret & 

273 Steinberg, 1977; Melrose, 1980; Bastian, 1994]. However, if the intrinsic source size and 9 C 

274 are similar in size, then the observed source size will be larger than 9 C by a factor near \/2. 

275 Note that this factor cannot account for the large discrepancies between the minimum 

276 source size observed at 327 MHz [Mercier et al., 2006] and those predicted here. If, on 
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277 the other hand, the intrinsic source size is much larger, then scatter broadening does not 

278 play a significant role. 

279 We have included refractive index effects that can be important when the radiation is 

280 emitted near the fundamental plasma level, but found them to be relatively unimportant. 

281 The inner scale is included via the Coles and Harmon [1986] model, interpreted in terms 

282 of cyclotron damping of MHD waves, and so depends primarily on the ambient electron 

283 density. We employ a hybrid model for the electron density that yields reasonable heights 

284 for meter wavelength emission at the fundamental. In view of the uncertainity in its value 

285 in the inner corona, the power law index a characterizing the turbulent spectrum is taken 

286 to be a free parameter. We consider both plane wave and spherical wave propagation. 

287 For the geometry we consider, where the source is embedded in the scattering medium, 

288 the spherical wave description is arguably more appropriate. 

289 We have thus explored a wide variety of effects. We observe that there is no significant 

290 difference in the predicted scattering angle between fundamental and second harmonic 

291 emission. We also find that the removal of the refractive index effect causes a negligible 

292 change in the predicted scattering angle. We find that the spherical divergence effect 

293 results in a significant lowering of the predicted scattering angle (by around 2 orders of 

294 magnitude). We find that removing inner scale effects by artificially setting the inner 

295 scale to be equal to a very small value (instead of determining it self-consistently from 

296 Eq 4) results in a significant enhancement of the predicted source size. The enhancement 

297 is greatest for flatter spectra, where it can be a factor of around 50, and it progressively 

298 disappears for steeper spectra. 
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99 As mentioned earlier, the power law index of the turbulent spectrum is a free parameter, 

ao There is a formal divergence at a = 4 in Eqs (6) and (7), and we therefore limit the 

01 computations to a maximum value of a = 3.97. The maximum value of the predicted 

02 scattering angle thus occurs at a — 3.97. 

03 For plane wave propagation, the predicted source size for Kolmogorov turbulence is 

04 around 10" lower than the observed one. For spherical wave propagation, we find that 

05 the maximum value of the predicted scattering angle is at least 25 times smaller than 

06 the observed one. It is emphasized that plane wave propagation is relevant to the well- 
o? accepted empirical formula of Erickson [1964], which predicts large source sizes, since it 
os pertains to observations of celestial sources through the solar wind. Even so, with current 

09 estimates of h(R) implying that s c ff > h (Figure 4), additional scattering is required to 

10 bring Erickson's result into quantitative agreement with the calculations here for plane 
a wave scattering. Alternatively, the inner scale k(R) should be smaller than that predicted 

12 by the Coles and Harmon [1989] model assumed here, so that s e g < U. 

13 The crucial result of this paper is that the predicted source sizes are considerably smaller 
« than the observed lower limit of 49" when spherical wave propagation and inner scale 
is effects are included, as they should be for sources in the solar corona. This broad trend of 
i6 the models substantially underpredicting the source size can be interpreted in three ways 
I? that are not exclusive and can occur in combination. First, it could imply that source 
is sizes much smaller than those that have been observed so far actually exist in the solar 

19 corona, and can potentially be observed. All the instances of observations of small sources 

20 to date have been limited by the instrument resolution; it is therefore quite likely that 

21 smaller sources can be detected when instrument resolutions are improved. Second, this 
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322 broad trend can be taken to imply that our naive extrapolation of the empirical form for 

323 the turbulence amplitude C%(R) to the inner corona is not justified. The results could 

324 be taken to imply that Cjj(R) in the inner corona is far higher than suggested by the 

325 empirical formula (3). This could be due to the functional form for C%(R) increasing 

326 more rapidly with decreasing R or due to a larger normalization factor or both effects. 

327 Third, the model (4) may significantly overestimate k(R), meaning that the turbulent 

328 cascade extends to smaller length scales (larger q) and leads to more scattering. These 

329 proposals all appear reasonable, and we regard all three as viable. 

330 Finally, we discuss the connection between our work and that of Bastian [1994]. The 

331 methodology is similar, and we investigate similar issues such as the effects of the inner 

332 scale, turbulence index a, and spherical versus planar wave propagation. Bastian's [1994] 

333 findings are contrary to ours: we find that our model predictions are substantially below 

334 the minimum observed size of 49", while Bastian's [1994] model predictions are substan- 

335 tially above 49". Thus, a priori, both models need revision. A major difference is in the 

336 choice of a model for C%. Bastian [1994] uses a model for C% which is proportional to 

337 the square of the background electron density and assumes the density model of Riddle 

338 (1974), which also involves a constant of proportionality. These two constants of propor- 

339 tionality are absorbed into one and fixed by normalizing the structure function D 2 o(10km) 

340 for a baseline of 10 km, an observing wavelength of 20 cm, and an elongation of 5 R & . In 

341 contrast, as explained earlier, the C% model we use is determined by an empirical fit to 

342 VLBI scattering observations between 10-50 R Q ; this was motivated by the need to use 

343 a C% model that is derived as directly as possible from observations. A minor matter is 

344 that Bastian [1994] discusses the disk to limb variation in the predicted scattering angle, 
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345 whereas our treatment is valid only for sources that are reasonably close to disk center. 

346 In order to do so, we would need to use the general formalism used here, together with an 

347 integration path that incorporates the appropriate extra path length needed for sources 

348 that are displaced from the disk center. 

349 An appropriate means of comparing the normalizations of the two treatments is thus 

350 to compare the normalizations of the structure function. Using (10) - (13), we write the 

351 structure functions as 



352 where 9 cs f(s) is the value of 9 C for spherical wave propagation using branch (7) and 

353 cp f(s) corresponds to plane wave propagation for branch (6). Then using our models for 

354 C%(R), k(R), and n hyb (R) we find that D sf (s) = 2.8 x 10~ 3 rad 2 for s = 10 km, A = 91 cm 

355 (corresponding to 327 MHz), a = 11/3 and a starting height corresponding to 327 MHz 

356 fundamental emission. The same prescription and parameters yield -D p /(s) = 7.7 x 10~ 3 

357 rad 5 / 3 . Since the structure functions we derive are based on integrations over heliocentric 

358 distance, we cannot assign a specific elongation to them. 

359 In comparison, Bastian normalizes by assuming D 2 ocm(10km) =4-12 rad 2 , based 

360 on measurements of by Coles and Harmon [1989] and Armstrong et al. [1990] of cosmic 

361 sources (implying primarily planar wave effects) at an elongation of 5 R Q . In order to 

362 normalize Bastian's [1994] values for the structure function to a wavelength of 91 cm, 

363 we concentrate on the structure function for spherical wave propagation. Inspection of 



364 Eqs [11] and [18] reveals that D sf (s) oc A 2 . Therefore, Z} 2 ocm(10km)/£> 91cm (10km) = 21. 

365 Bastian's [1994] range of values for D 2 ocm(10km) thus corresponds to -D 9 i cm (10km) = 




(18) 
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see 82 — 250 rad 2 . The difference in D 9 i cm (10km) between the two prescriptions is thus a 

367 factor of (30 — 90) x 10 3 , corresponding to a factor m 170 — 300 in 9 C . 

sea This large difference ~ (30 — 90) x 10 3 in the normalization of the structure function 

369 is primarily indicative of a corresponding difference in the normalization of C% between 

370 the two treatments; this is because neglect of inner scale effects increases 6 C by less than 

371 a factor of 10 for Kolmogorov turbulence in Figure 7. In this connection, we note that 

372 Bastian's [1994] normalization for D 2 ocm(10km) is based on values of D(s) measured at 

373 an elongation of 5 R®. We also note (e.g., Fig 1 of Coles & Harmon [1989]) that values 

374 of D(s) measured at larger elongations can be considerably lower (by as much as a few 

375 orders or magnitude, depending upon the elongation). This is significant, since the model 

376 for Clf that we use in this paper is based on observations between 10 and 50 R Q . 

377 In summary, the foregoing results demonstrate conclusively that spherical wave prop- 

378 agation effects are vital for solar sources, with plane wave predictions several orders of 

379 magnitude larger than the spherical predictions. Similarly, inner scale effects are quan- 

380 titatively important, while fundamental versus harmonic radiation effects are relatively 

381 small. The results and discussion above demonstrate the importance of accurate models 

382 for Cjf(R) and to a lesser extent models of U{R) and a(R). This paper's prescription for 

383 Cj; (Eq [3]) is empirical and directly based on observations (but extrapolated to smaller 

384 R), does not have any normalization constants that need to be determined, and leads to 
sas scattered sizes for a point source that are smaller than the minimum source size observed 
386 to date (49" by Mercier et al. [2006]). Thus smaller source sizes than 49" may be observ- 
ss? able. In contrast, another well-known prescription [Bastian, 1994] predicts much stronger 
388 scattering with source sizes always larger than 49": while this is inconsistent with the 
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sag minimum source size observed to date at 327 MHz, it may provide the extra scattering 

390 required to account for Erickson's empirical angular broadening result for cosmic sources 

391 viewed through the solar wind. 

392 While this paper's results extend and confirm previous theoretical results pertaining to 

393 spherical vs. plane wave effects and provide the first explanation of the small source sizes 

394 recently observed, it is also clear that more observational and theoretical work is required 

395 on Cjf(R) especially, but also on k(R) and a(k, R). This includes temporal variations over 

396 the solar cycle but also spatial variations between radio source regions and other regions 

397 of the corona. Increases in Cfj(R) and decreases in U(R) would lead to more scattering. 

398 Work on both n e (R) and 5n(R)/n e (R) may be useful [e.g., Efimov et al., 2008; Cairns et 

399 al, 2009]. It is quite possible that scattering observations and theory will provide useful 
wo constraints on these five quantities and therefore on the processes heating the solar corona 

401 and accelerating the solar wind. 
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Figure 1. Geometry for spherically diverging wavefront, which is appropriate to a 
situation where the source is embedded in the scattering medium 
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Figure 2. The plasma frequencies predicted by the four density models in the text 
are plotted against the height h = R — R Q above the photosphere (in units of Rq). The 
solid line uses the Leblanc et al model (15), the dotted line uses the Aschwanden & Benz 
model (16), the dashed line uses the 4*Newkirk model (14) and the dash-dot line uses the 
"hybrid" model (17). 
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Figure 3. The inner scales (in units of km) predicted by some density models in the 
text are plotted against the height h = R — R Q above the photosphere (in units of R Q ). 
The solid line uses the Leblanc et al model (15), the dashed line uses the 4*Newkirk model 
(14), and the dash-dot line uses the hybrid model (17). 
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Figure 4. The ratio s e ffA is plotted against the height h = R — R Q above the 
photosphere (in units of it!©). The solid line shows the ratio for plane wave propagation 
(i.e., with s c ff = s), while the dotted line shows the ratio for spherical wave propagation 
(i.e., with s e ff = sR/{R\ — Rq)). The hybrid density profile is used for determining the 
density and the inner scale. 
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Figure 5. Predicted 9 C (1 AU) in at 327 MHz for plane wave propagation, as a 
function of a, using (6), (10) and (13). The observed source size is 49". The thin line is 
for fundamental emission and the thick line is for second harmonic emission. 
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Figure 6. Predicted # c (1 AU) in " at 327 MHz for plane wave propagation as a function 
of a, using Eqs (7), (10) and (13). The solid line is for fundamental emission and the 
dotted line is for second harmonic emission. The heavy lines are computed with Zj set at 
an artificially low value of 1 m, while still using branch (7). 



DRAFT 



February 8, 2011, 4:13am 



DRAFT 



X - 34 SUBRAMANIAN & CAIRNS: CONSTRAINING CORONAL TURBULENCE WITH RADIO DATA 



N 

*(0 

(D 

O 
\- 

o 

(f) 

TD 

0) 

O 



"O 

0) 



10' 



10 ( 



10 



-2 



^ 10 



10 



-4 



-6 




3.0 



Spherical wave 



S e ff << "i 



3.2 



3.4 3.6 
a 



a = 11/3 
Kolmogorov 



3.8 



4.0 



Figure 7. Predicted # c (1 AU) in " at 327 MHz for spherical wave propagation as 
a function of a, using (7), (11) and (13). The observed source size at 327 MHz is 49". 
There is negligible difference between fundamental and second harmonic emission. The 
solid lines are computed with U from prescription (4), while the dotted lines are computed 
with li set to an artificially low value of 1 m, while still using branch (7). 
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Figure 8. Predicted # c (1 AU) in " at 327 MHz for spherical wave propagation as 
a function of a, using (6), (11) and (13).. The difference between the predictions for 
fundamental and second harmonic emission is negligible. 
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Figure 9. Direct comparison of 9 C (1 AU) for spherical (solid lines) and plane wave 
(dotted lines) propagation. Fundamental emision is assumed. 
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